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Phot on- phot on Absorption of Very High Energy Gamma-Rays 
from Microquasars: Application to LS 5039 

Markus Bottcher 1 & Charles D. Dermer 2 
ABSTRACT 

Very high energy (VHE) 7-rays have recently been detected from the Galactic 
black-hole candidate and microquasar LS 5039. A plausible site for the produc- 
tion of these VHE 7-rays is the region close to the base of the mildly relativistic 
outflow. However, at distances comparable to the binary separation, the intense 
photon field of the stellar companion leads to substantial 77 absorption of VHE 
7-rays. If the system is viewed at a substantial inclination (i ^ 0), this absorp- 
tion feature will be modulated on the orbital period of the binary as a result 
of a phase-dependent stellar-radiation intensity and pair-production threshold. 
We apply our results to LS 5039 and find that (1) 77 absorption effects will be 
substantial if the photon production site is located at a distance from the central 
compact object of the order of the binary separation (pa 2.5 x 10 12 cm) or less; 
(2) the 77 absorption depth will be largest at a few hundred GeV, leading to a 
characteristic absorption trough; (3) the 77 absorption feature will be strongly 
modulated on the orbital period, characterized by a spectral hardening accom- 
panying periodic dips of the VHE 7-ray flux; and (4) 7-rays can escape virtually 
unabsorbed, even from within pa 10 12 cm, when the star is located behind the 
production site as seen by the observer. 
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1. Introduction 

Recent observations (Aharonian et al. 2005) of > 250 GeV 7-rays with the High Energy 
Stereoscopic System (HESS) from the X-ray binary jet source LS 5039 establish that micro- 
quasars are a new class of 7-ray emitting sources. These results confirm the earlier tentative 
identification of LS 5039 with the EGRET source 3EG J1824-1514 (Paredes et al. 2000). In 
addition to LS 5039, the high-mass X-ray binary LSI 61°303 (V615 Cas) also has a possible 
7-ray counterpart in the MeV - GeV energy range (Gregory & Taylor 1978; Taylor et al. 
1992; Kniffen et al. 1997). 

Micro quasars now join blazar AGNs as a firmly established class of very-high energy 
( VHE; with energies > a few hundred GeV - TeV) 7-ray sources. The nonthermal continuum 
emission of blazars is believed to be produced in a relativistic plasma jet oriented at a small 
angle with respect to our line of sight. Their radio through UV/X-ray emission is most likely 
due to synchrotron emission by relativistic electrons in the jet, while the high energy emission 
can be produced by Compton upscattering of lower-energy photons off relativistic electrons 
(for a recent review, see, e.g. Bottcher 2002), or through hadronic processes (Mannheim & 
Biermann 1992; Atoyan & Dermer 2001; Miicke et al. 2003). 

Because of their apparent similarity with their supermassive AGN cousins, it has been 
suggested that Galactic microquasars may be promising sites of VHE 7-ray production (e.g., 
Bosch- Ramon et al. 2005a). While earlier work on 7-ray emission from X-ray binaries focused 
on neutron star magnetospheres as 7-ray production sites (e.g. Moskalenko et al. 1993; 
Moskalenko & Karakula 1994; Bednarek 1997, 2000), the VHE 7-ray detection of LS 5039 
suggests that the 7-ray production is more likely to be associated with the jets. High-energy 
7-rays of microquasars can be produced via hadronic (e.g. Romero et al. 2003) or leptonic 
processes. In the latter case, the most plausible site would be close to the base of the 
mildly relativistic jets, where ultrarelativistic electrons can Compton upscatter soft photons. 
Possible sources of soft photons are the synchrotron radiation produced in the jet by the 
same ultrarelativistic electron population (SSC = synchrotron self-Compton; Aharonian & 
Atoyan 1999), or external photon fields (Bosch- Ramon & Paredes 2004; Bosch- Ramon et al. 
2005a). Both LS 5039 and V615 Cas are high-mass X-ray binaries which are rather faint in 
X-rays, with characteristic 1-10 keV luminosities of ~ 10 34 ergs s -1 . This is much lower 
than the characteristic bolometric luminosity of the high-mass companions of these objects, 
at > 10 38 erg/s. Consequently, the dominant source of external photons in LS 5039 and 
V615 Cas is the companion's optical/UV photon field. 

The intense radiation field of the high-mass companion will, however, also lead to 77 
absorption of VHE 7-rays in the ~ 100 GeV - TeV photon energy range if VHE photons 
are produced close to the base of the jet. For the case of LS 5039, Aharonian et al. (2005) 
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have estimated that VHE emission produced within ~ 10 12 cm of the central engine of this 
microquasar would be heavily attenuated (r 77 ~ 20 for E 1 ~ 100 GeV), and the observed 
spectrum of VHE photons would be hardened compared to its intrinsic shape. For com- 
parison, the 77 opacity due to stellar radiation back-scattered by the stellar wind in the 
vicinity of the binary system (for which the angle of incidence between the line of sight and 
the target soft photon field would be more favorable for 77 absorption) may be estimated as 
r™*{r) ~ (L* Mcrf,)/(2.7 • 48 7r 2 r 2 m p m e c 3 e*), where is the stellar luminosity, M is 
the mass outflow rate, is the terminal velocity of the stellar wind, and e* is the dimension- 
less mean photon energy of the stellar radiation field. For LS 5039, ~ 7 x 10 38 ergs s -1 , 
m e c 2 e* ~ 3.5 eV, M « 10- 6 - 3 M©/yr, and v<x « 2500 km s" 1 (McSwain & Gies 2002), 
which yields r" 7 ax (r) ~ 0.1/r 2 2 , where r = 10 12 ri2 cm. Consequently, the 77 opacity will 
be strongly dominated by the direct stellar photon field. In the same spirit, one can also 
estimate the effect of 7-ray absorption in the field of atomic nuclei in the stellar wind. Using 
a total absorption cross section per unit mass of <jh ~ 0.012 cm 2 g _1 for GeV - TeV 7-rays, 
we find r 7 z PS er# M/(47ri>oo r ) ~ 1.3 x 10~ 4 / r i2 for LS 5093, which is also always much 
smaller than the 77 absorption depth. 

In this Letter, we provide a more detailed analysis of the expected 77 absorption trough 
caused by direct companion star light, including its temporal modulation due to the orbital 
motion. The model description and the derivation of the general expression for the 77 
opacity as a function of 7-ray photon energy and orbital phase is given in §2. Numerical 
results for LS 5039 are presented in § 3. § 4 contains a brief summary and discussion of our 
results. 



2. Model description and analysis 

We choose a generic model set-up as illustrated in Fig. 1 (for a comparable model setup 
in the AGN case see Bottcher & Dermer 1995). The orbital plane of the binary system 
defines the (x±, x 2 ) plane. The jet, assumed to be perpendicular to this plane, defines the x 3 
axis. The system is inclined with respect to our line of sight by an inclination angle i. An 
azimuthal (phase) angle is defined such that <fi = in the direction of the X2 axis. The line 
of sight lies in the (x 2 , £3) plane. The 7-ray production site is located at a height zq along 
the jet. The phase angle of the companion star at the time of production of a 7-ray photon 
at z is denoted 0o- The distance the photon travels along the line of sight is denoted by 
I, whereas a starlight photon travels a distance x before interacting with the 7-ray photon. 
The angle of incidence between the two photons is 9, and /j, = cos 9. The star was located at 
azimuthal angle <f>i at the time when a photon, leaving the star at that time, interacts with a 
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Fig. 1. — Geometry of the model. The direction of the radio jets defines the £3 axis. The 
orbital plane of the binary system is the (xi, x 2 ) plane, defined in such a way that line of 
sight towards the observer lies in the (x 2 , £3) plane, where the azimuthal angle = 0. 

7-ray photon that has travelled a distance / from zq. This yields 0i = 0o + (27r/P) (I — x)/ c, 
where P is the orbital period. 

With the definitions as shown in Fig. 1, the distance x can be calculated as x = |x| = 
|r — s|, where s is the orbital separation of the binary system. We find 

x 2 = s 2 + I 2 + Zq + 2 / (z cosi — s sini cos0i). (1) 

From Fig. 1, /x = 1 • x/(/x), which yields 



/ + Zq cos i — s sin i cos 
H = 



x 



(2) 
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With these quantities, we can evaluate the 77 opacity of a 7-ray photon with dimen- 
sionless energy e = E /(m e c 2 ) as 

00 00 

r 77 (e , zo, <j> ) = J dl (1 - fi) J de(T^(e ,e,fi)n* ph (e,x) (3) 
2 

where <r 77 is the pair production cross section and 

n P h(e ' X) = 4n 5 m e c 3 Qix 2 (e £ / e *-l)' (4) 

The stellar spectrum is approximated by a blackbody with dimensionless temperature O* = 
kT*/m e c 2 , and the star is approximated as a point source. 



3. Results for LS 5039 

Casares et al. (2005) have analyzed optical intermediate-dispersion spectroscopic ob- 
servations of LS 5039, which provided an estimate of the mass of the compact object of 
Mx = 3.7lJ;oM Q , indicating that it is likely to be a black hole. Other relevant parameters 
of the binary system are P = 3.91 d, the luminosity of the 06.5V type stellar companion, 
= 1O 5 - 3 L , with an effective surface temperature of T cS = 39,000°K, a viewing angle of 
i = 25°, and an orbital separation of s ~ 2.5 x 10 12 cm (Casares et al. 2005). Note that 
the center of mass is rather close to the massive (M rs 23M q ) stellar companion. However, 
because the light-travel and jet propagation times on the length scales investigated here 
(< 10 14 cm) are much shorter than the orbital period, the detailed dynamics of the orbital 
motion are irrelevant for our purposes: at any given phase <f>o, the geometry is basically 
stationary. An additional complication could be introduced by the substantial eccentricity 
of the orbit, e = 0.35, with periastron at a phase angle of ~ 0.6 tt (Casares et al. 2005). 
However, in view of the still rather poor quality of the data currently available and expected 
in the near future, this may be regarded as a higher-order effect which may be introduced 
at a later stage. 

The inset to Fig. 2 illustrates the shape of the absorption trough caused by 77 absorption 
and its dependence on the orbital phase. Here we assume that the intrinsic 7-ray spectrum is 
a power-law with photon index a ph = 2.5, and z = 10 12 cm. The various curves illustrate the 
orbital modulation of the absorption trough, with the lowest (most heavily absorbed) curve 
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corresponding to O = and the highest (least absorbed) curve corresponding to (fi = n. 
The modulation is a combined consequence of two effects: for phase angles closer to n, (a) 
the average distance of the star to any point on the line of sight is longer and (b) the angle 
of incidence 9 is smaller. In addition to effect (a) causing the overall photon number density 
of the stellar photon field to decrease, effect (b) causes the threshold for 77 pair production 
to increase as e t hr = 2/(e* [1 — //]). This leads to a decreasing overall depth of the absorption 
trough, and a shift of the minimum of the absorption trough towards higher photon energies. 

Fig. 2 shows the dependence of the absorption feature on the location zq of the VHE 
7-ray production site. The 77 opacity is plotted for two photon energies, E = 250 GeV, 
and E — 1 TeV at O — 0- GeV - TeV photons produced within z ~ s from the compact 
object will be heavily attenuated for this phase angle. For photons produced at z 3> s, 
77 attenuation becomes negligible. The crossing-over of the two curves in Fig. 2 illustrates 
the shift of the minimum of the absorption trough to higher photon energies with increasing 
distance from the central engine. This is a consequence of the decreasing incidence angle as 
discussed above. 

Fig. 3 quantifies the orbital modulation of the 77 absorption trough for the same pa- 
rameters as used in the inset to Fig. 2. The dotted and dot-dashed curves show the or- 
bital modulation of the optical depth, illustrating that (a) it is maximized for 0o ~ 0, 27r, 
and (b) that the mimimum of the absorption trough shifts towards higher photon ener- 
gies for phases closer to ir. The most feasible way to detect an orbital modulation would 
be to measure a periodicity of the observed photon flux. The solid and short-dashed 
curves in Fig. 3 show the integrated photon number flux above energies Eq = 250 GeV 
and Eq — 1 TeV, respectively. This illustrates that the relative modulation is most pro- 
nounced at several hundred GeV, noting that the 77 absorption is largest at photon energies 
Ei = 2 • 2(m e c 2 ) 2 / (kT^l — cos(arctan(s/,2 ) +*)]). In this estimate, we have made use of the 
fact that the photon energy t pk at the peak of the 77 attenuation cross section is e pk 2e thr . 
For the case of LS 5039 with zq = 10 12 cm and 0o — 0, E\ = 300 GeV, in good agreement 
with the inset in Fig. 2. The spectral softening with increasing flux is also illustrated by the 
local photon spectral index 0:500 at 500 GeV, plotted as the long-dashed curve in Fig. 3. The 
flux dip around 0o ~ is accompanied by a hardening of the VHE 7-ray spectrum. 

Because the bulk speed of microquasar jets is typically only mildly relativistic (r ~ 2), 
in dramatic contrast to blazars, one might also need to consider the possibility of VHE 7- 
rays emanating from the counter jet. Thus, we have done a study as described above, for an 
inclination angle of % = 155°. While the results for this case are qualitatively very similar 
to the i = 25° case, the resulting 77 opacities are typically larger by a factor of ~ 2 - 3 for 
photon production sites at z < 10 12 cm, as expected because of the more favorable photon 
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collision angle and the line of sight passing by the star closer than for the approaching-jet 
case. Because of the stronger absorption of 7 rays from the counter jet throughout the binary 
orbit, and the reduction in the flux due to Doppler deboosting, the counter jet is unlikely to 
enhance the 77 absorption signature in the signal appreciably. As for % = 25°, 77 absorption 
effects become negligible for zo ^> s. 

4. Summary and Discussion 

We have presented a detailed analysis of the effect of 77 absorption of VHE 7-rays 
near the base of the jet of a microquasar by the intense photon field of a high-mass stellar 
companion. We include the time- dependent, periodic modulation of this effect due to the 
binary's orbital motion. We applied our results to the specific case of LS 5039, which has 
recently been identified as the counterpart of the VHE 7-ray source HESS J1826-148. Our 
results can be summarized as follows: 

(1) VHE 7-rays produced closer to the central engine than zq ~ a few xlO 12 cm, which 
is of the order of the binary separation s, would be subject to very strong 77 absorption due 
to the stellar radiation field at orbital phases close to 0o — 0. 

(2) For VHE photon production sites at z < s, the 77 opacity - - and, thus, the 
VHE 7-ray flux — would be strongly modulated on the orbital period of the binary system 
(P = 3.91 d in the case of LS 5039). At orbital phases close to 0o = tt, the intrinsic VHE 
7-ray flux would still be virtually unabsorbed even for Zo ~ 10 12 cm. 

(3) The orbital modulation of the VHE 7-ray flux would be characterized by a spectral 
hardening in the ~ 300 GeV - 1 TeV range during flux dips. At lower energies, the spectrum 
softens with decreasing flux. 

The HESS collaboration has reported no evidence for either flaring or periodic variability 
from LS 5039 (Aharonian et al. 2005), although Casares et al. (2005) suggest that there 
is weak evidence for variability of the HESS emission with the orbital period. Periodic 
variability may also be indicated by X-ray observations (Bosch-Ramon et al. 2005b) of LS 
5039. 

Besides 77 opacity effects, there are a few alternative scenarios which might cause a 
periodic modulation of the 7-ray flux: 

(a) As a consequence of the substantial eccentricity of the orbit, the rate of mass transfer 
from the stellar companion to the compact object, which is believed to be dominated by wind 
accretion, is likely to be periodically modulated. This modulation would also be expected 
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to appear at radio and X-ray energies. 

(b) Analogous to the phase-dependent modulation of the incidence angle for 77 absorp- 
tion, this geometrical effect would also yield a more favorable angle for Compton scattering 
of starlight photons into the 7-ray regime at phases 0o ~ 0. 

(c) The orientation of the jet may also be mis-aligned with respect to the normal of 
the orbital plane (Maccarone 2002; Butt 2003) and possibly precessing about the normal 
(Larwood 1998; Torres et al. 2005), leading to additional modulations, including a changing 
Doppler boosting factor. 

Effect (a) would be expected to lead to an overall hardening of the 7-ray spectrum at all 
energies with increasing 7-ray flux, while (b) would lead to an overall softening throughout 
the GeV - TeV photon energy range because of the Klein-Nishina cutoff becoming noticeable 
at lower observed photon energies with increasing flux. Both effects are in contrast to the 
77 absorption trough investigated in this Letter. A stationary misalignment of the jet could 
lead to a slight enhancement of the orbital modulation (if the jet makes a smaller angle with 
the line of sight than the orbital- motion axis) or reduce it (in the opposite case), but would 
not change our results qualitatively. A 7-ray flux modulation due to jet precession can easily 
be disentangled from the orbital modulation since the precession period is generally different 
from the orbital period, so that its effect would average out when folding observational data 
with the orbital period. Consequently, a measurement of a non-thermal absorption trough at 
VHE 7-ray energies modulated with the orbital period would firmly establish the importance 
of 77 absorption effects and thus place a robust limit on the distance z of the VHE 7-ray 
production site in LS 5039. 
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Fig. 2. — 77 opacity at 250 GeV and 1 TeV as a function of the distance of the photon 
production region from the central compact object at phase <po = 0. The figure illustrates 
that (1) VHE photons produced within a few xlO 12 cm (i.e., of the order of the orbital 
separation of the binary system) would be subject to substantial 77 absorption; (2) the 
minimum of the absorption trough (maximum of function of photon energy) is 

shifting towards higher energies for larger distances from the central source. Inset: Orbital 
modulation of the expected 77 absorption trough, assuming a power-law spectrum with 
photon index a ph = 2.5 and a photon production site at z = 10 12 cm. The different curves 
represent the escaping photon spectrum at various orbital phases, from O — (lowest curve) 
to 0o = Ti" (highest curve) in steps of n/10. 
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Fig. 3. — Orbital modulation of the integrated photon number flux above energies E = 
250 GeV (solid) and E = 1 TeV (short-dashed), the 77 opacity at E = 250 GeV (dotted) 
and E = 1 TeV (dot-dashed), and the local photon spectral index 0:500 at 500 GeV (long- 
dashed). As in the inset to Fig. 2, an underlying power-law of photon index o p h = 2.5 and 
a photon production site at z = 10 12 cm has been assumed. A periodic flux modulation is 
expected to be accompanied by positive spectral-index/flux correlation (spectral softening 
as the flux increases) at E > 300 GeV; the opposite behavior is found at E < 100 GeV. 



